Stability and thermodynamics of brane black holes 
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I. INTRODUCTION 

The extra dimensions idea had its origin in the seminal 
works by Kaluza and Klein in the 20's dj and gained mo- 
mentum in the context of string theory in the last decades 
[3- Recent developments on higher-dimensional grav- 
ity resulted in a number of interesting theoretical ideas 
such as the brane world concept. The essence of this 
string inspired model is that Standard Model fields are 
confined to a three dimensional hypersurface, the brane, 
while gravity propagates in the full spacetime, the bulk. 
The simplest models in this context (abbreviated here as 
RSI and RSII), proposed by Randall and Sundrum Q, 
describe our world as a domain wall embedded in a Z2- 
symmetric five dimensional anti-de Sitter (AdS) space- 
time. The RSI model proposes a mechanism to solve the 
hierarchy problem by a small extra dimension, while the 
RSII model considers an infinite extra dimension with a 
warp factor which ensures the localization of gravity on 
our brane. 

Black holes are important sources of gravitational 
waves that are expected to be detected by the cur- 
rent and upcoming generation of experiments. This will 
open up a new window for testing modifications of gen- 
eral relativity. The simplest case of gravitational col- 
lapse in the standard four dimensional world is described 
by the 4-dimensional Schwarzschild metric. In the 5- 
dimensional scenario it would be natural to ask whether 
matter confined on the brane after undergoing gravita- 
tional collapse can still be described by a Schwarzschild- 
type metric. The most natural generalization in the 
RSII model corresponds to a black string infinite in the 
fifth dimension, whose induced metric on the brane is 
purely Schwarzschild 4]. However, although the cur- 



vature scalars are everywhere finite, the Kretschmann 
scalar diverges at the AdS horizon at infinity, which turns 
the above solution into a physically unsuitable object. 

It has been argued that there exists a localized black 
cigar solution with a finite extension along the extra di- 
mension due to a Gregory-Lafiamme 5] type of instabil- 
ity near the AdS horizon. A class of such a solution has 
been found by Casadio et al. 0, Q using the projected 
Einstein equations on the brane derived by Shiromizu et 
al. 8]. It has the desired "pancake" horizon structure 
ensuring a non-singular behavior in the curvature and 
Kretschmann scalars at least until the order of the mul- 
tipole expansion considered there. In fact, this solution 
belongs to a class of black hole solutions found later by 
Bronnikov et al. |3|, who also classified several possi- 
ble brane black holes obtained from the Shiromizu et al. 
projected equations 8J in two families according to the 
horizon order. For such spacetimes only horizons of order 
1 or 2 are possible, but not higher than that. 

In this paper we are interested in the study of black 
holes from the point of view of a brane observer, as our- 
selves. We analyze some characteristics of Bronnikov et 
al. solutions. Firstly, some aspects of the thermody- 
namics are studied. Black Hole Thermodynamics was 
constructed when Bekenstein proposed the proportional- 
ity law between the entropy and the horizon area |inl |. 
The discovery of Hawking radiation validated this pro- 
posal and established the proportionality factor 1/4 in 
a precise way |llj leading to the well-known Bekenstein- 
Hawking formula. 
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One way to compute the entropy based on a semi- 
classical description of a scalar field was proposed by 't 
Hooft ,12] and it is known as the brick wall method, which 
was frequently used later in several contexts [I^l . When 
applying this method to a Schwarzschild black hole, 't 
Hooft found that the entropy was proportional to the 
area, as expected, but additionally it had a a^^ correc- 
tion, a being the proper distance from the horizon to 



the wall. This term was later interpreted as a one-loop 
correction to the Bekenstein-Hawking formula, since it 
can be absorbed as a renormalization of the gravitational 
couphng constant G [l^. 

Another interesting feature of black hole thermody- 
namics is the existence of an upper bound on the entropy 
of any neutral system of energy E and maximal radius 
R in the form S < 2TrER, proposed by Bekenstein [l5l |. 
This bound becomes necessary in order to enforce the 
generalized second law of thermodynamics (GSL). 

Besides thermodynamical results, we also consider the 
response of a brane black hole perturbation which should 
represent some damped oscillating signal. It can be de- 
composed with Laplace transformation techniques into a 
set of so-called quasinormal modes (QNM). The QNMs 
of black holes are important because they dominate in 
the intermediate late-time decay of a perturbation and 
do not depend upon the way they were excited. They 
depend only on the parameters of a black hole and are, 
therefore, the "footprints" of this structure. The time- 
independent problem for perturbations of a brane black 
hole turns out to be quite similar to that for a black hole: 
one has to find the solutions of the wave-like equations 
satisfying the appropriate boundary conditions, which we 
shall further discuss in detail. 

The paper is organized as follows, in Sec. II the brane 
black holes are presented. Sec. Ill discusses the ther- 
modynamical properties of the solutions thus considered. 
Sec. IV treats the question of perturbation and stability 
of these objects. Wave-like equations for the perturba- 
tions are derived. In Sec. V an analysis of the quasinor- 
mal modes and late-time behavior is developed. In Sec. 
VI we summarize our results, and some final comments 
are made. 



II. BRANE BLACK HOLE SOLUTIONS 

The vacuum Einstein equations in 5 dimensions, when 
projected on a 4-dimensional spacetime and after intro- 
ducing gaussian normal coordinates (a;'', with /n = . . . 3, 
and z) , lead to the gravitational equation on the 3-brane 
given by Q 



R^n} = A45lt^ - E^. 



(2) 



where A4 is the brane cosmological constant, and -B^^ is 
proportional to the (traceless) projection on the brane of 
the 5-dimensional Weyl tensor. 

The only combination of the Einstein equations in a 
brane world that can be written unambiguously without 



specifying E^i, is their trace [a, Q > 

i?(4) = 4A4 



(3) 



It is clear that this equation, also known as the hamilto- 
nian constraint in the ADM decomposition of the met- 
ric, is a weaker restriction than the purely 4-dimensional 



equation Rfj,^, = 0, which, in fact, is equivalent to Eq. 
(O, provided that we know the structure of E^^^. 

In order to obtain four dimensional solutions of Eq. 
(O, we choose the spherically symmetric form of the 4- 
dimensional metric given by 



ds' 



-A{r)dt^ 






r^ide^ +iim'^ 01 



(4) 



We relax the condition A{r) = B{r), which is acciden- 
tally verified in four dimensions but, in fact, there is no 
reason for it to continue to be valid in this scenario. In 
this spirit, black hole and wormhole solutions 6, 9, 16, ij 
as well as star solutions on the brane [13 have been ob- 
tained in the last years. We should mention that even 
without relaxing this condition, previous solutions have 
also been found JI3, |23 ■ 

In this context the hamiltonian constraint can be writ- 
ten explicitly in terms of A and B as 



B 



A^ 
~A 



' 1 A! B' 

+ 2AB+' 

;{l-B)=AAi 



B^ 



A^ 
~A 



(5) 



with prime (') denoting differentiation with respect to r. 
We will center our attention in the black hole type 
solutions which can be obtained by one of the following 
algorithms BHl and BH2. These are subclasses of the 
corresponding algorithms in [3 (where the parameter s 
in this reference is set to 1). 

First Algorithm (BHl): Specify a function A{r), positive 
and analytical in a neighborhood of the event horizon 
Tl[r], in such a way that 4 A + rA' > in 7?.[r], and 
A fa (r — r/i) as r — > r/j. Then B{r) is given by the 
general solution of ((SJ with vanishing brane cosmological 
constant, 



Bir) = 



2Ae^^ 



r{4:A + rA')^ 
X f {4:A + rA'){2-r^R)e-^^dr + C 

-Jrh 

where C is an integration constant and 



r(r) 



A' 



4A 4- rA! 



■ dr . 



(6) 



(7) 



For C > we have a black hole metric with a horizon at 
r = rji, which is simple if C > and of the order 2 + p if 
both C = 0, and Q has the behavior. 



Qir) 



r^R: 



{r — rhY , near r = r/j , p € N . (8) 



Second Algorithm (BH2): Specify a function A(r), posi- 
tive and analytical in a neighborhood of the event hori- 
zon 7?.[r], in such a way that A A + rA' > in 7?.[r], and 
A K (r — rfi)'^ as r — > r/i. Then B{r) is again given by 
(0. The black hole metric appears when C = with 



a horizon a,t r = rh of the order 2 + p ii Q{r) behaves 
according to l|HJl- 

Both algorithms lead to double horizons in the case 
C = 0, iiQ{rh) > 0. A case in point of the first algorithm 
is the solution with the metric element A{r) having the 
usual form of a Schwarzschild black hole found by Casa- 
dio, Fabbri, and Mazzacurati (CFM solution) [a, 13 given 

by 



Mr) 
B{r) 



1 



2M 



(l-2M)(l-i^) 



1 _ 3M 
^ It 



(9) 



where 7 is an integration constant. The event horizon is 
localized at r^ = 2M and the singularity at r = 2>M/2 
instead of r = 0. Notice that the Schwarzschild solution 
is recovered with 7 = 3. In this work we are restricted to 
the case when 7 < 4. 

Another interesting example of this algorithm is the 
metric with zero Schwarzschild mass given by 



A{r) 
B{r) 



1- — , h>0, 






1 



C-h 

V2r2 - h^ 



(10) 



whose horizon r = ft, is simple if C > and double if 
C = 0. The singularity occurs at r = h/y/2. This exam- 
ple shows that in the brane world context a black hole 
may exist without matter and without mass, only as a 
tidal effect from the bulk gravity. However, there is a spe- 
cial situation when h^ can be related to a 5-dimensional 
mass, namely, C — h. In this case Eq. (|10|l is the induced 
metric of a 5-dimensional Schwarzschild black hole, as de- 
scribed in ,21,], where the chosen background was ADD- 
type. 



III. BLACK HOLE THERMODYNAMICS 

In order to study the thermodynamical properties of 
the brane black holes generated by the BHl and BII2 
algorithms, we use the following expressions of the metric 
coefficients near the horizon 

A{r) = Ai{r-rh)+0{{r-rHf) (11) 

B{r) - B^C{r-rh) + B2{r-rhf + 0{{r-rhf), 

for BHl algorithm with Ai,Bi,B2 > and C being an 
integration constant that defines the black hole family, 
and 

A{r) = A2{r - rh)^ + 0{r - rh)\ 

B{r) = B^C + B^{r - Thf + B^{r - rhf 

+0{{r-ruf), (12) 



We will show here the calculation for the BH2 family, 
which turns out to be more interesting, since the metric 
coefficients expansion (|12|l is different from the standard 
one (|11|) . However, we will display certain quantities for 
both families wherever it is relevant. 

We first consider the issue of the entropy bound. The 
surface gravity at the event horizon is given by 



iV^iBiC for BHl family, 
y/A^B^ for BH2 family . 



(13) 



Let us consider an object with rest mass m and proper 
radius R descending into a BH2 black hola The con- 
stants of motion associated to t and are |22J 



E = TTt , J = — 7r0 , 



where 



In addition, 



TTt 



9ttt, 



(14) 



(15) 



(16) 



For simplicity we just consider the equatorial motion of 
the object, i.e., 9 =^ ix 12. The quadratic equation for the 
conserved energy E of the body coming from 114|l - ll()|) is 
given by 



aE^ - 2PE + C = 0, 



with 



a 

c 



(17) 



(18) 



A2ir~rhy{j' + m'r^) 



The gradual approach to the black hole must stop when 
the proper distance from the body's center of mass to the 
black hole horizon equals R, the body's radius. 



Th+SiR) 



dr 



R. 



(19) 



VW) 

Integrating this equation we obtain the expression for S, 



j \CR^Bi for BHl family , 
I Ry/B^ for BH2 family . 



(20) 



Solving H17|l for the energy and evaluating at the point 
of capture r = Vh + 5 we have 



Er. 



^A2{.P+m'^rl)5 
rh 



(21) 



This energy is minimal for a minimal increase in the black 
hole surface area, J — 0, such that 



for BH2 algorithm, being C the family parameter again. 



E„ 



' Ai m 5 . 



(22) 



From the First Law of Black Hole Thermodynamics we 
know that 



Now we introduce the semiclassical quantization con- 
dition, 



dM = -dAr , 
2 ' 



(23) 



where A^ is the rationalized area (Area/47r), and dM = 
Emin is the change in the black hole mass due to the 
assimilation of the body. Thus, using H13I) we obtain 



dAr = 2mi? . 



(24) 



Assuming the validity of the GSL, Sbh{M + dM) > 
Sbh{M) + S*, we derive an upper bound to the entropy 
S of an arbitrary system of proper energy E, 



S < 2TrER . 



(25) 



This result coincides with that obtained for the purely 
4-dimensional Schwarzschild solution, and it is also inde- 
pendent of the black hole parameters [l^ • It shows that 
the bulk does not affect the universality of the entropy 
hound. 

Let us find now the quantum corrections to the classi- 
cal BH entropy. We consider a massive scalar field $ in 
the background of a BH2 black hole satisfying the mas- 
sive Klein-Gordon equation, 



□ - 



^)$==0. 



(26) 



In order to quantize this scalar field we adopt the sta- 
tistical mechanical approach using the partition function 
Z, whose leading contribution comes from the classical 
solutions of the euclidean lagrangian that leads to the 
Bekenstein-Hawking formula. In order to compute the 
quantum corrections due to the scalar field we use the 't 
Hooft's brick wall method, which introduces an ultravio- 
let cutoff near the horizon, such that 

$(r) = at r = rh+£, (27) 

and an infrared cutoff very far away from the horizon, 

$(r) = at r = L-»M. (28) 

Thus, using the black hole metric Q) and the Ansatz 
~ R{r)Ye,nid,<l>), F,q.^ turns out to be 



(J, ^ e^'^* 



E^ 

A 



R 



Wf>^ 


{r' 


[ r2 ^ 


m^ 



\/ABdrR\ 
R = 0. 



(29) 



Using a first order WKB approximation with R{r) w 
giS(r) jj^ IJ29JI and taking the real part of this equation we 
can obtain the radial wave number K = drS as being. 



i^ = i?-V2 



E^ 
A 



fi{i+ll 
\ r2 



1/2 



(30) 



K{r,£,E)dr. 



(31) 



rh+e 



In order to compute the entropy of the system we first 
calculate the Helmholtz free energy F of a thermal bath 
of scalar particles with temperature 1/P, 

F=^ f d£{2l + 1) / dur ln(l - e"''^) . (32) 

Integrating by parts, using Ip^Hjl and lp?T|l. and performing 
the integral in £ we have 



37r 



F = -A / drA-^/^B-i/V dE 



rh+e 



{E^ - Am2)3/2 



rJ3E_l 



(33) 

Following 't Hooft's method, the contribution of this 
integral near the horizon is given by 



2r' 



,3 /•£ 



dy^^^^'^y"^-'^-' r dE- 



E^ 



^~ 37r .A^."" (i?.3a)i/2 7, -^e/3^-1' 

(34) 

where y = r/rh, s = e/rh, and L = L/rh- Notice that 
since A goes to near the horizon, the mass term in Eq. 
(|33|l becomes negligible. 

Therefore, the leading divergent contribution to F 
(with £ ^ 0) is 



r^3 (^^)-3/2 

The corresponding entropy is then, 

^dF _ 4r2 7r3(A2)-3/2 



Se = f3' 



d(3 45(B3C)i/2£2^3 



(35) 



(36) 



Using the value of the Hawking temperature Th = 1//3 = 
k/2tt, with K given in (|13|l we have 



5,= 



rlBiC /{360e) for BHl family , 

rlBsC / {90 e"^) for BH2 family . ^ ' 



We can express our result in terms of the proper thick- 
ness a given by 

+" dr ^ J2^/y/B^ for BHl family , 
^/B{r) ~ [e/VW^ for BH2 family . 

(38) 



Thus, 



S,= 



90a2 



or in terms of the horizon area Area = 47rr2, 



S, 



Area 
360^2 



(39) 



(40) 



which is the same quadratically divergent correction 
found by 't Hooft .121 for the Schwarzschild black hole 
and by Nandi et al. jl^ for the CFM brane black hole. 
Thus, we see that the correction is linearly dependent on 
the area. 

The calculation of the entropy bound and entropy 
quantum correction for the BHl black hole is similar and 
leads to the same results shown in 1^5)1 and ^^. 



IV. PERTURBATIVE DYNAMICS: MATTER 
AND GRAVITATIONAL PERTURBATIONS 

For simplicity we model the matter field by a scalar 
field $ confined on the brane obeying the massless 
{m = 0) version of the Klein-Gordon equation H26(l . We 
expect that massive fields (m ^ Oj should show rather 
different tail behavior, but such cases will not be treated 
in the present paper. 

Using the decomposition of the scalar field as 
'i'{t,r,9,(j)) = i?(i,r)Yf_m(0, 0) in terms of the angular, 
radial, and time variables we have the equation 



d^Ri 



+ ■ 



VsMr.))Ri, 



dt'^ dri 

with the tortoise coordinate r^ defined as 
dr^,{r) 1 



(41) 



(42) 



dr y/A{r)B{r) 

The effective potential Vgc is given by 

T4, = A(r)fci^ + -^ [Air) B' [r) + A'{r) B{r)] . 
j,z 2r 

(43) 
In order to address the problem of black hole stability 
under gravitational perturbations, we consider first order 
perturbation of R^p = —Eap, where R^ff and Eap are 
the Ricci tensor and the projection of the five dimensional 
Weyl tensor on the brane, respectively. In general, the 
gravitational perturbations depend on the tidal pertur- 
bations, namely, SEa/s- Since the complete bulk solution 
is not known, we shall use the simplifying assumption 
5Eap = 0. This assumption can be justified at least in 
a regime where the perturbation energy does not exceed 
the threshold of the Kaluza-Klein massive modes. Anal- 
ysis of gravitational shortcuts [23| also supports this sim- 
plification showing that gravitational fields do not travel 
deep into the bulk. On the other hand, since we ignore 
bulk back-reaction, the developed perturbative analysis 
should not describe the late-time behavior of gravita- 
tional perturbations. Within such premises we obtain 
the gravitational perturbation equation 



6Ra.B = . 



(44) 



We will consider axial perturbations in the brane ge- 
ometry, following the treatment in (23. To make this 
section more self contained, we will briefly describe the 



treatment used in [23. The metric Ansatz with sufficient 
generality is 



ds^ 



2 - e"'dt'-e^''(^^^^ 



— Lodt — q2d9 — qadcf) 



-e'^'^^dr^ - e'^'^^de'^ (45) 

where we adopt here a more convenient notation, 
1 



e^-^ - A{r) , e^^^ = -^ , e^^^ = r^ , e^^^ = r^ sin^ 9 . 
Bir) 

(46) 
Axial perturbations in the brane metric I0J are charac- 
terized by non-null (but first order) values for tu, q2, and 
qs in Eq. l|15|) . We refer to [23 for further details. 

In order to decouple the system, it is adopted the 
change of variables 



and 



Qap — qa,l3 — qi3,a , 



QaO ~ Qa,0 ^ t^iQ ) 



(47) 



(48) 



with a,f3 = 2,3. We denote partial differentiation with 
respect to t, 9 and by ",0", ",1" and ",2", respectively. 
The perturbations are then described by 

= -(^,2 - g2,o),o , (49) 



Cg3l/' + !^-A'2-A»3Q2Q) „e~3'/' + !^ + A«3-A»2 _ 

= (^^,3 - <73,o),o • (50) 

Setting (5(i,r, 6*) = cxp{3tjj + iy- fi2- t^3)Q23, Eqs.lgni) 
and H5U|I can be combined as 



lB{r) d 
A{r) dr 



^^A{r)B{r) 



sm 



dQ_ 
dr 

1 dQ 
89 \sm'^9 89 



lB{r)8^Q 

AW at2" 



(51) 



We further separate variables and write Ea. (|51() in the 
form of a Schrodinger-type wave equation by introduc- 
ing Q{t,r,9) =rZi{t,r)C^j^^ (9), and r = r{ri,), where 

^e+2 (^) ^^ ^^^'^ Gegenbauer function. Thus, the axial 
gravitational perturbations are given by an equation of 
motion with the form given in 141|l with the effective po- 
tential 



V.r..ir)-Air)^I±^^^'-^^ 
-l.[A{r)B'ir)+A'ir)Bir)]. 



(52) 



SPECIFIC MODELS 



Overview of the results 



The equations of motion for the scalar and axial grav- 
itational perturbations give us a tool to analyze the dy- 
namics and stability of the black hole solutions in both 
the CFM and "zero mass" black hole backgrounds. Of 
particular interest are the quasinormal modes. They are 
defined as the solutions of Eq. 141() which satisfy both 
boundary conditions that require purely out-going waves 
at (brane) spatial infinity and purely in-going waves at 
the event horizon, 



lim *e=^'"^ 



constant , 



(53) 



with ^ = i?£ and Z^ for the scalar and gravitational 
perturbations, respectively. 

In order to analyze quasinormal mode phase and late- 
time behavior of the perturbations, we apply a numerical 
characteristic integration scheme based in the light-cone 
variables u = t — r^, and v = t + r^, used, for exam- 
ple, in [2^ [23, |23|. In addition, to check some results 
obtained in "time-dependent" approach we employ the 
semi-analytical WKB-type method developed in 29] and 
improved in |33| • Both approaches show good agreement 
for the fundamental overtone which is the dominating 
contribution in the signal for intermediate late-time. 

At a qualitative level we have observed the usual pic- 
ture in the perturbativc dynamics for all fields and ge- 
ometries considered here. After the initial transient 
regime, it follows the quasinormal mode phase and finally 
a power-law tail. In contrast to the 5-dimensional model 
in |31j . in the present context we do not observe Kaluza- 
Klein massive modes in the late-time behavior of the per- 
turbations. This is actually expected, since our treat- 
ment for gravitational perturbations neglects the back- 
reaction from the bulk. Still, as discussed in section IV, 
our results should model the quasinormal regime. 

A necessary condition for the stability of the geome- 
tries we have considered is that these backgrounds must 
be stable under the perturbations modelled by the effec- 
tive potentials l|IS|) and ifS^ . 

If the effective potential (V) is positive definite, the 
differential operator 
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Figure 1: (Left) Effective potential for the scalar perturba- 
tions in the CFM background for very negative values of 7. 
Negative peaks are displayed in detail. (Right) Bounded evo- 
lution of the scalar field perturbation with such effective po- 
tentials. The dotted line is the late-time power-law tail. The 
parameters are £ — and AI — 1. 
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Figure 2: (Left) Effective potential for the axial gravitational 
perturbations in the CFM background for very negative values 
of 7. (Right) Bounded evolution of the gravitational field 
perturbation with such effective potentials. The dotted line 
is the late-time power-law tail. The parameters in the graphs 
are ^ = 2 and M = 1. 



V 



(54) 



is a positive self-adjoint operator in the Hilbert space 
of square integrable functions of r^,, and, therefore, all 
solutions of the perturbativc equations of motion with 
compact support initial conditions are bounded. 

However, as we will see, the effective potentials may 
be non-positive definite for certain choices of the param- 
eters in Eas.(j!!|l- HlU|) . Nevertheless, even when the effec- 
tive potential is not positive definite, we do not observe 
unbounded solutions. 



Using both high order WKB method and direct nu- 
merical integration of the equations of motion a numeri- 
cal search for quasinormal modes with positive imaginary 
part was performed for scalar and gravitational pertur- 
bations. One of the most important results in this work 
is that no unstable mode was observed. Furthermore, the 
perturbativc late-time tails have power-law behavior (in 
one case an oscillatory decay with power-law envelope) . 



Table I: Fundamental quasinormal frequencies for the scalar perturbation around the CFM black hole for several values of the 
parameters 7 and £. The black hole mass is set to M = 1. 
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Re(a;o) 


-lm(ajo) 


Re(ajo) 


-Im(c<Jo) 


Re(c^o) 


-Im(c<Jo) 




-5 


0.28580 


0.19779 


0.208204 


0.225080 


0.305499 


0.181027 




-2 


0.29201 


0.16608 


0.276143 


0.181776 


0.309350 


0.122589 







0.29337 


0.14138 


0.299359 


0.154590 


0.236213 


0.165871 




1 


0.29400 


0.12853 


0.298384 


0.134468 


0.252449 


0.171879 




2 


0.29415 


0.11378 


0.294679 


0.115208 


0.293168 


0.120632 




3 


0.29283 


0.098045 


0.291114 


0.0980014 


0.292910 


0.0977616 




3.9 


0.29076 


0.082598 


0.287181 


0.0820285 


0.289628 


0.0811812 


2 


-5 


0.48053 


0.18951 


0.488726 


0.207518 


0.433498 


0.183580 


2 


-2 


0.48266 


0.16069 


0.488518 


0.166137 


0.451478 


0.193079 


2 





0.48413 


0.13815 


0.486195 


0.139258 


0.485925 


0.142572 


2 


1 


0.48449 


0.12570 


0.485289 


0.126043 


0.485911 


0.126758 


2 


2 


0.48447 


0.11206 


0.484420 


0.112159 


0.484691 


0.112283 


2 


3 


0.48317 


0.097097 


0.483211 


0.09680485 


0.483642 


0.0967661 


2 


3.9 


0.48178 


0.080778 


0.481091 


0.08098300 


0.481705 


0.0808983 



Table II: Fundamental quasinormal frequencies for the axial gravitational perturbation around the CFM black hole for several 
values of the parameters 7 and i. The black hole mass is set to M = 1 . 







Direct Int 


sgration 


WKB-3"' 


order 


WKB-e**^ 


order 


£ 


7 


Re(cjo) 


-Im(ci;o) 


Re{uJo) 


-Im(a;o) 


Ro(c^o) 


-Im(c<Jo) 


2 


-5 


0.36409 


0.18017 


0.401345 


0.197274 


0.418575 


0.193276 


2 


-2 


0.37049 


0.16062 


0.391442 


0.163223 


0.402747 


0.163222 


2 





0.37359 


0.13539 


0.384611 


0.137147 


0.389781 


0.139044 


2 


1 


0.37457 


0.12138 


0.381053 


0.122624 


0.383017 


0.124656 


2 


2 


0.37483 


0.10604 


0.377306 


0.106767 


0.377126 


0.107996 


2 


3 


0.37368 


0.088957 


0.373162 


0.0892174 


0.373619 


0.0888910 


2 


3.9 


0.36961 


0.072435 


0.368552 


0.0717786 


0.371935 


0.0712303 


3 


-5 


0.59476 


0.18365 


0.608026 


0.191340 


0.613069 


0.194241 


3 


-2 


0.59835 


0.15845 


0.604567 


0.159482 


0.605032 


0.161682 


3 





0.59993 


0.13527 


0.602594 


0.135600 


0.601901 


0.136695 


3 


1 


0.60033 


0.12238 


0.601646 


0.122520 


0.601033 


0.123070 


3 


2 


0.60026 


0.10832 


0.600614 


0.108374 


0.600375 


0.108525 


3 


3 


0.59947 


0.092690 


0.599265 


0.0927284 


0.599443 


0.0927025 


3 


3.9 


0.59700 


0.077176 


0.597227 


0.0767434 


0.597584 


0.0767411 



B. CFM brane black holes 



We first consider scalar perturbations in the CFM sce- 
nario. The tortoise coordinate r* after the explicit inte- 
gration is 



T2(r) = ^^^^ ^' In [4r - M(3 + 7) + 2Ti(r)] , (57) 



r,{r)^Ti{r) + T2{r)+T3{r) 



(55) 



n{r)=- 



2M 



with 



In 



M(5 - 7)r - M^(6 - 7) -h My/4^Ti{r) 
r-2M 



Ti (r) = - v/(2r-7M)(2r-3M) , (56) 



(58) 
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The scalar and axial gravitational effective potentials 
for perturbations in the CFM background (respectively, 
V^^^ and V^^^) in terms of the parameters M and 7 
are given by 



M(7- 



2M' 

1 

r 

3)(r2 -6Mr 



1) 2M 



6Af 2) 



r3(2r-3Af)2 



(59) 



and 



yCFM 



(r) 



2M 
r 



6M 



Af (7 - 3)(5r2 - 20Afr + ISAf^) 



r3(2r-3A/)2 



(60) 



Setting 7 = 3 we recover the usual expressions for pertur- 
bations around the four dimensional Schwarzschild black 
hole. 

A basic feature of the effective potentials V^^^' and 
^grav'^ is that they are not positive definite. Typically, 
for negative enough values of the parameter 7 (with large 
I7I) a negative peak in the effective potential will show 
up. It is no longer obvious that the scalar and gravita- 
tional perturbations will be stable. 

If the effective potential is not positive definite (and 
cannot be approximated by a positive definite one), the 
WKB semi-analytical formulae will usually not be appli- 
cable, but direct integration techniques will. Using the 
later approach an extensive search for unstable solutions 
was made. One important result in this work is that even 
for very high values of 7, the (scalar and gravitational) 
perturbative dynamics is always stable. This is illustrated 
in Figs. nand[21 where we show non-positive definite ef- 
fective potentials and the corresponding (bounded) field 
evolution. 

The overall picture of the perturbative dynamics for 
the effective potentials hereby considered is the usual 
one. After a brief transient regime, the quasinormal 
mode dominated phase follows, and, finally, at late times 
a power-law tail dominates. 

For the fundamental multipole mode {£ — 0) the effec- 
tive potential will not be positive definite for any value 
of the parameter 7. Direct integration shows that the 
field evolution is always bounded for a great range of 
variation of 7. This point is illustrated in Figs. ^ and 
2, where we have selected rather large values of 7. In- 
deed, it is observed that the decay is dominated from very 
early time by the power-law tail. Therefore, it is very 
difficult to estimate the quasinormal frequencies directly 
from this "time-dependent" approach. The WKB-type 
expressions are not applicable if £ = for two reasons: 
the effective potentials are not positive for r larger than 
a certain value, and it is well known that this method 
works better with £ < n, where n is the overtone num- 
ber. 

With small but non-zero values of £ the quasinormal 
frequencies can be accurately estimated. As it is shown in 



Tables I and II, the concordance with the WKB results is 
reasonable, except for some values of 7 (typically around 
7 = 0). 

For large values of £ an analytical expression for the 
quasinormal frequencies can be obtained. Expanding the 
effective potential in terms of small values of l/£ and 
using the WKB method in the lowest order (which is 
exact in this limit) we find 



Re(w„) 



3V3M 



Im(cj„ 




(61) 



(62) 



As it can be seen from the data in tables I and II 
and Fig. |31 the dependence of the frequencies with the 
parameter 7 is very weak, although not trivial. In a 
large variation range of 7 the absolute value of Im(wo) is 
a monotonically decreasing function, while Re(wo) typi- 
cally has maximum points. 

The late-time behavior of the perturbations considered 
here can be treated analytically. Far from the black hole 
the scalar effective potential in terms of r^ assumes the 
form, 



2M 



v: 



CFM 



(r.) 



-1) , 4Me(e+l)ln{rj2M} 

I 1 ;3 



with £=^0 

with £>0. 

__ " (63) 

It is then shown [3j, |33| that with the initial data having 
compact support a potential with this form has a late- 
time tail 



r: 



CFM 



t' 



-(2^+3) 



(64) 



Therefore, at asymptotically late times the perturbation 
decays as a power-law tail for any value of the parameter 
7. This is a strong indication that the models are indeed 
stable. This point is illustrated in Figs, ^andl^l It is 
reminiscent from a similar behavior of the Gauss Bon- 
net term added to Einstein gravity in higher dimensions, 
which was recently treated in |3J|. Although the result 
is formally valid also for the gravitational perturbations 
we have considered, it should be noted that in the sim- 
plified model developed in this paper the back-reaction 
from the bulk, which can modify the tail presented here, 
was neglected. 



C. "Zero mass" brane black holes 

The treatment in section IV is general enough to in- 
clude also the case of perturbations around the "zero 
mass" brane black hole. Using the metric (|10|l the 
scalar and axial gravitational perturbations are de- 
scribed by wave equations similar to Ea. H41|l with 



Table III: Fundamental quasinormal frequencies for the scalar perturbation around the "zero mass" black hole for several values 
of the C and I. The parameter h is set to h = 1. 







Direct Integ 


Tation 


WKB-3"' 


order 


WKB-6*'' 


order 


I 


C 


Re(a;o) 


-Im(a;o) 


Re(a;o) 


-Im([iJo) 


Re(a;o) 


-lm(tJo) 











0.728358 


0.232800 


0.746504 


0.230697 




0.1 


0.74494 


0.24551 


0.730716 


0.247607 


0.749727 


0.245155 




0.5 


0.75196 


0.30135 


0.732795 


0.303801 


0.755273 


0.301191 




0.9 


0.75437 


0.35023 


0.726005 


0.356124 


0.751242 


0.357060 




1.0 


0.75410 


0.36201 


0.722899 


0.368853 


0.748461 


0.371463 




1.1 


0.75354 


0.37352 


0.719182 


0.381471 


0.745143 


0.385894 




2.0 


0.74548 


0.46711 


0.660519 


0.494129 


0.732352 


0.479133 


2 





_ 


_ 


1.242071 


0.230667 


1.246220 


0.231061 


2 


0.1 


1.24534 


0.21001 


1.243863 


0.245695 


1.248248 


0.246014 


2 


0.5 


1.25534 


0.28301 


1.246290 


0.299897 


1.251377 


0.300033 


2 


0.9 


1.25072 


0.34681 


1.244715 


0.347044 


1.250452 


0.347062 


2 


1.0 


1.24989 


0.35677 


1.243914 


0.358055 


1.249839 


0.358066 


2 


1.1 


1.25377 


0.36031 


1.242968 


0.368816 


1.249090 


0.368835 


2 


2.0 


1.24523 


0.44843 


1.227403 


0.456844 


1.235779 


0.459396 



Pi 





n = 




0.615 






\ 1 1 1 1 1 




0.610 


-\ 


- 


0.605 


\^ 


- 


0.600 





N 



-10-8 -6-4-2 2 4 





and 



-10-8 -6-4-2 2 4 

Y 



-10-8 -6-4-2 2 4 

Y 



Figure 3: Dependence of gravitational perturbation quasinor- 
mal frequencies on 7 in the CFM geometry. The results are 
qualitatively similar for the scalar perturbation. The param- 
eters are i = 2> and M — 1. 



effective potentials given by 



vr'MHi-pf'''"" ''' 



C-h 



(2r2 - /l2)3/2 



1- 



r^ r 



(65) 



^ grav 



+ 



(r) 



_h^ 



(2r2 - /i2)3/2 



5- 



£{£ + !) 

11 /i2 4h' 

1 

.4 



4^ 

r4 






r* 



(66) 



Again, the effective potentials can be non-positive defi- 
nite for specific choices of parameters, as illustrated in 
Fig. H (left). For example, Hi = and C > h, F/™ will 
not be positive definite. If ^ > 0, VfJ^ and Vg^^v ^^^^ be 
non-positive definite for high enough values of C. 

Except for C — h, an explicit expression for the tor- 
toise coordinate was not found. Nevertheless, the nu- 
merical integration is possible. The semi-analytical WKB 
approach was also used to compute quasinormal frequen- 
cies. The concordance is excellent. The WKB formulas 
seem to be more reliable in the present case. With the 
choice C = 1 we recover some results considered in [21[ . 
Again, an extensive search for unstable modes was per- 
formed. Some calculated frequencies are shown in Tables 
III and IV. Our results show that the dynamics of the 
scalar and axial gravitational perturbations is always sta- 
ble in the "zero mass" background. We illustrate this 
point in Figs. ^and|5l 

Analytical expressions for the quasinormal frequencies 
for the scalar and gravitational perturbations can be ob- 
tained in the limit of large multipole index £. As done in 
the CFM geometry, we obtain 



Re(w„) 



2/1 



(67) 
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Figure 4: (Left) Effective potential for the scalar perturba- 
tions in the "zero mass" black hole background with high 
values of C . Negative peaks are displayed in detail. (Right) 
Bounded evolution of the scalar field perturbation with such 
effective potentials. The dotted line is the late-time power-law 
tail. The parameters are i = and h = 1. 
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Figure 5: (Left) Typical effective potential for the gravi- 
tational perturbations in the "zero mass" black hole back- 
ground. (Right) Bounded evolution of the gravitational field 
perturbation with such effective potentials. The parameters 
are 1 = 2 and h = 1. 



Ini(w„) 



V2/1 



(68) 



We expect that the late-time behavior of the gravita- 
tional perturbations should be dominated by back scat- 
tering from the bulk, not considered here. But the 
tail contribution to the scalar decay can be analytically 
treated, at least in the limit where r ^ h. In this case 
ii C > h or < C < h, the effective potential T^c™ is 
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Figure 6: Bounded evolution of the scalar field perturbation in 
the "zero mass" black hole background with C = h for several 
values of £. After the quasinormal mode phase a power-law 
tail is observed. The tail dependence with £ obeys Eq. 11731 . 
In the graphs the parameter h was set to ^ = 1. 



approximated by 



v^rir*) 



■ C-h 
V2r3 



with £ = , 

_ _ (69) 

Again, we observe that ([3j, yj]) with the initial data 
having compact support the tail has the form 



R'r - t-(2f+3) ^j^jj 0<C <h or Oh. 



(70) 



An interesting limit is when C = h. In this case the 
explicit expression for the tortoise coordinate according 
to the usual definition in Eq. (|42|l is 



h (r \ h 
r*(r) =r + -ln - - 1 - -In 



1 



(71) 



The effective potential V^^™ with C — h is approximated 
by 



v^rir.) 




2l{l+l)h^ 



with £ = , 
with £>0. 



(72) 



In this limit a power-law still dominates the late-time 
decay. But its dependence with the multipole index £ is 
different, 



Rf 



f-(2m) ^itl^ c^h. 



(73) 



This point is illustrated in Fig. 

As observed in the CFM model, for the non-extreme 
"zero mass" model the scalar perturbation decays as a 
power-law tail suggesting that the model is stable. 

The qualitative picture of the field evolution in the 
"zero mass" black hole — quasinormal mode followed by 
power-law tail — changes drastically when the extreme 
case (C — 0) is considered (see Fig. 0). If £ = 0, we ob- 
serve the usual power-law tail dominating the late-time 
decay. But when £ > 0, the simple power-law tail is re- 
placed by an oscillatory decay with a power-law envelope, 



Rf 



t~^/^ sin {ui X t) with C = and £ > . (74) 



11 











10 


10 






















,i 


1 




ur 


! 


\ 




v^ 




10"* 


- 


\ 




1^ 


itei' ' 


IIV" 








"in 


k 




! 


I 


= 3 


" 'W 


ii 


.„-s 


: 




ii 


' 


™ 



Figure 7: Bounded evolution of the scalar field perturbation 
in the extreme "zero mass" black hole background (C = 0), 
for several values of £. If ^ = 0, the decay is dominated by 
a power- law tail (i~^). If ^ > 0, the decay is dominated by 
an oscillatory tail, whose envelope is t~ ' . In the graphs the 
parameter h was set to h — 1. 
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Figure 8: (Left) Numerical value of t where the scalar 
wave function is zero in the extreme "zero mass" black hole. 
Straight lines imply that the period of oscillation is a con- 
stant. The parameter h is set to ft = 1. (Right) Dependence 
of the angular frequency uji in Eg. 1741 with I for several values 
of ft. 



Therefore, for £ > the power index (—3/2) is inde- 
pendent of the multipole index (.. The angular frequency 
LOi for large times approaches a constant, as we can see 
in Fig. |Sl(left). The angular frequency is well approxi- 
mated by a linear function of £, as indicated in Table V for 
some values of h. This result implies that the dominating 
contributions in the late-time decay are the modes with 
C. > Q, i.e., the power-law enveloped oscillatory terms. 
We also observe that these tails dominate from very early 
times, so that it was not possible to estimate the quasi- 
normal frequencies in the extreme case (as indicated in 
the first lines of Tables III and IV). 



VI. CONCLUSIONS 

In this work we studied brane black holes from the 
point of view of a brane observer. We considered the 
two family solutions found by Bronnikov et al. |9| in 
order to derive the Bekenstein entropy bound and the 
one-loop correction to the Bekenstein-Hawking formula 
using the 't Hooft brickwall method. In addition, we 
performed scalar and axial gravitational perturbations 
in two specific examples of these families. With these 
perturbations we were able to analyze the dynamics and 
stability of the black hole solutions. 

The results of the black hole thermodynamics study 
show that the entropy bound continues to be independent 
of the black hole parameters. Thus, the presence of the 
bulk does not affect the universality of the entropy bound 
for a brane observer, as ourselves, reinforcing the Gen- 
eralized Second Law. Moreover, applying the 't Hooft's 
brickwall method to both black hole families we see that 
the entropy correction takes the same form as that of 
a Schwarzschild black hole when written in terms of its 
own black hole parameters. Therefore, as the correction 
is linearly dependent on the area, it can be absorbed in 
a renormalized gravitational constant. 

One of the most important results in this paper came 
from the perturbative dynamics. We should stress that 
the assumption SEa/3 = was necessary in order to solve 
the gravitational perturbation equation (|44|) without any 
knowledge of the bulk structure. This vanishing tidal ef- 
fect is perfectly justified when the perturbation energy 
is lower than the threshold of the Kaluza-Klein massive 
modes. Likewise, as we neglect the bulk back-reaction, 
our analysis does not describe the perturbation late-time 
behavior. Our results show no unstable mode in the 
scalar and gravitational analysis. In addition, the late- 
time tails display a power-law behavior what enforces 
their stability. 

In the case of CFM black hole even if the effective 
potential is not positive-definite the quasinormal modes 
are stable (negative imaginary part). The agreement of 
the several methods employed in the calculation is good 
for i not too small. 

On the other hand, in the case of the "zero mass" black 
hole we observe a richer picture. The scalar and gravita- 
tional field evolution is always bounded suggesting that 
this class of models is stable. But the late time decay 
of the matter field strongly depends on the parameters 
C and h. If C is non-zero and not equal to /i, the late- 
time decay is dominated by a power-law tail with the 
usual dependence on the multipole parameter i. But if 
C — h, this dependence changes. Finally, in the extreme 
regime (C = 0) the late-time decay is dominated by os- 
cillatory modes with a power-law envelope. This power 
index seems to be universal, not depending on £. 

Summarizing, the thermodynamics in the class 
of models we considered is consistent, while 
the dynamics in specific backgrounds is stable 
in the approach employed in the present work. 
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Table IV: Fundamental quasinormal frequencies for the axial gravitational perturbation around the "zero mass" black hole for 
several values of the C and I. The parameter h is set to /i = 1. 







Direct Integ 


;ration 


WKB-3"' 


order 


WKB-6*'' 


order 


I 


C 


Ro(c^o) 


-Im(a;o) 


Re{LOo) 


-Im(c<Jo) 


Re(a;o) 


-Im(a;o) 


2 









0.934530 


0.191198 


0.934386 


0.219631 


2 


0.1 


0.94412 


0.21924 


0.935381 


0.205749 


0.964066 


0.215448 


2 


0.5 


0.94783 


0.26533 


0.938066 


0.264005 


0.958654 


0.249590 


2 


0.9 


0.94334 


0.30949 


0.942035 


0.318153 


0.929971 


0.317851 


2 


1.0 


0.94207 


0.31993 


0.943255 


0.330875 


0.928449 


0.334460 


2 


1.1 


0.94069 


0.33016 


0.944551 


0.343278 


0.928937 


0.350243 


2 


2.0 


0.92442 


0.41343 


0.958486 


0.442675 


0.973998 


0.456121 


3 





_ 


_ 


1.537001 


0.214671 


1.538619 


0.215204 


3 


0.1 


1.50122 


0.25435 


1.537512 


0.229197 


1.539283 


0.229786 


3 


0.5 


1.53958 


0.27762 


1.536138 


0.280829 


1.538697 


0.280900 


3 


0.9 


1.53905 


0.32341 


1.532735 


0.325500 


1.534597 


0.324529 


3 


1.0 


1.53757 


0.33402 


1.531809 


0.335949 


1.533070 


0.334820 


3 


1.1 


1.53594 


0.34411 


1.530883 


0.346171 


1.531369 


0.344973 


3 


2.0 


1.52161 


0.42100 


1.523732 


0.430806 


1.513596 


0.433952 



Table V: Oscillatory frequency of the tail in the extreme 
mass" (C = 0) black hole for several values of h. 



h 



Angular Frequency ijJi {I > 0) 



0.5 
1.0 
2.0 



0.02593 + 0.06421 x t 
0.01347 + 0.03191 x^ 
0.006683 + 0.01597 x£ 



While our results suggest that the brane models 
presented are viable, the final check would be the 



analysis of the continuation in the bulk of the geometries 
presented here. 
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